STRUCTURE OF BINARY SEQUENCES 



J. THARRATS 

Abstract. The distribution of a given sequence in the set of all sequences with n ones and 
777 = M — n zeros are found by relating the problem to the partitions of a natural number 
in m natural summands, taking into account the order. The formulas obtained have many 
applications both in Physics and Mathematics. Examples discussed in the present paper are: 
non Markovian chains, partition functions of binary alloys and Ising magnets, generalized 
Kaplansky lemma, generalized Fibonacci numbers and a general expansion of XX=o ^ (™) 
in terms of the Stirling numbers of second kind. 



1. Introduction 

Some combinatorial problems can be solved immediately by refering them to the structure 
of the ( m ^ n ) sequences of m zeros and n ones of length N = m+n, i.e.: knowing the ocurrence 
distribution of a given ordered string of r O's and s l's (r + s < N). Take, for example the 
problem of finding the number of subsets of cardinality n of the set X = {1, 2, ... , N} that 
contain exactly h times a number r of consecutive naturals; this is done knowing how many 
times the string (1, 1, . r . ., 1) has h recurrences in those sequences. For a non fixed cardinality 
one makes the sum over n and obtains how many elements of the power set have the said 
property. In this case, as one can see in Riordan || problem 1(b), p. 14], for h = 0, r = 2 
one obtains the Fibonacci numbers, though Carlitz || had generalized to h = 0, r > 2; (in 
1976 Carlitz [|[ precises the problem mixing strings of (00) and (11)). 

Thus, the structure of digit-sequences can have its own place in Combinatorics, introducing 
problems concerning the sequences themselves; for instance, in the (^) sequences, the number 
of sequences that have each pair of ones separated by at least k zeros. Although the solution 
is given by Kaplansky lemma || , the problem becomes trivial by applying the formula which 
gives the distribution of the (0, 0, 0; 1) strings, choosing the result which gives exactly n 
strings. More generally, we will also find the case in which only s ones will be separated by 
at least k zeros. 

In this theory two points of view can be taken: One may decide to bound the sequences 
with first and last digit or to consider the first digit following the last in cyclic fashion. Since 
our main interests are physical applications, only the second point of view will be followed. 
This view, called here "Ising process" gives in the particular case above the corrected Fi- 
bonacci numbers. It is also more symmetrical, which explains why our method is so simple 
and can be developed free of generating functions, in contrast with the method of bounded 
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sequences in which the particular cases solved by Carlitz give results strongly dependent on 
the first digit. 

This paper was motivated by the solution of non-markovian random walks.Q As is well- 
known, the solution of the one-dimensional random walk, with probabilities p (moving to 
the right) and q = 1 — p (to the left), is given by Bernoulli's distribution: 




But if transition probabilities are dependent on past states, i.e., if the particle has memory, 
then the distribution changes radically. In the simple case of dependence on the last state 
only, the probability is ap or aq for particle moving backwards and f3p or f3q moving forwards 
(a + j3 = 1). Let T™ n be the number of paths of length N with r changes in direction of 
the particle, or "jumps" in N steps (number of ocurrences of (01) or (10)). Then paths like 
ppqqp and pqpqp have different probability because they are affected by the factors a 2 (3 3 and 
a 4 /3 respectively, and the binomial coefficient in ([!]) splits in the form 




Thus the formula (|XJ) changes depending on the constant a which measures the tendency to 
go forwards (a < |) or backwards (a > |). We see that the problem is solved by finding the 
T-numbers. 

We also will see how T-numbers give the natural way to find some Thermodynamic Parti- 
tion Functions, such as in Ising magnets, in contrast with naive ways given by some physicists. 

Historically, it seems Andre Q was the first to consider such sequences (see Comptet f|) 
in order to solve Bertrand's scrutiny problem. He considered minimal paths in a square 
lattice where ?/-jumps stand for 0's and x-jumps for l's. This give a simple way to solve 
Bertrand's problem counting the number of paths which do not touch the y = x line. The 
solution is not invariant by a circular permutation of the sequences and falls in the case of 
unbounded sequences, for the path (0010001) is allowable (in each step the eventual winner 
is winning) but the Ising equivalent path (0100010) is not. Processes like this one are also 
non-markovian because allowed decisions depend of all anterior states, but can be solved by 
elementary methods; the scrutiny problem can be solved as in p] or by reducing it to well 
known random walks avoiding the origin (projecting paths on y = —x). 

Along the way in this paper many combinatorial formulas arise, some very well known 
but others related to Stirling numbers are new, as far as we know; but what really matters 
is the natural way in which they arise. 

Notation. We prefer to preserve the (old fashioned?) symbols for: C™ = ( n ) ; S™ = Stirling 
numbers of the second kind; P™ = number of partitions of a natural number n in m natural 
summands disregarding the order; M™ for number of partitions regarding the order (M 
stands for De Moivre who discovered it to be equal to ( n ~ 1 l )). 

1 The paper "Analytical Solution of non-Markovian Random Walks" in which we give some applications 
to crystal lattices, is to appear. 
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As P-partitions can be associated with Young-Sylvester's tableaux, M-partitions too, tak- 
ing into account the order. For example, the partition of 6 into 3 parts gives the Pg = 3 
tableaux: 

rrm rrp m 



For M-partitions, we have the following Mg = ( 2 ) = + jjfj u + §j = 10 tableaux: 



which we call Young-De Moivre tableaux. 

In the unbounded poset of integers, we call mAn the min (m, n) and mVn the max (m, n); 
(no confusion arise with connectors of formal systems). 

We will refer always the known combinatorial formulas to the monumental "Ars" by Knuth 

0-0 

2. T-NUMBERS 

Let us consider the set of all sequences with n ones and m = N — n zeros; an element of 
this set is, for example (for N = 24, m = 13): 

000001110011010001100111. (2) 

To every sequence we associate a number r, which is the number of jumps from to 1, or 
from 1 to 0, including extremes, so in the case of the sequence (0) we have r = 10. Obviously, 
the number r is even and it is invariant under any circular permutation of the sequence. 

The problem here is to find the distribution of r's, i.e.: How many r equal 2, how many 
t equal 4, etc. in these (^) sequences. 

The following table shows the case N = 7 and m = 3 with the corresponding values of r: 

Table I. 





T 




r 




r 




T 




T 


0000111 


2 


0011001 


4 


0101010 


6 


1000101 


4 


1010100 


6 


0001011 


4 


0011010 


4 


0101100 


4 


1000110 


4 


1011000 


4 


0001101 


4 


0011100 


2 


0110001 


4 


1001001 


4 


1100001 


2 


0001110 


2 


0100011 


4 


0110010 


4 


1001010 


6 


1100010 


4 


0010011 


4 


0100101 


6 


0110100 


4 


1001100 


4 


1100100 


4 


0010101 


6 


0100110 


4 


0111000 


2 


1010001 


4 


1101000 


4 


0010110 


4 


0101001 


6 


1000011 


2 


1010010 


6 


1110000 


2 



The table gives 7, 21 and 7 sequences (7 + 21 + 7 = Q) for the corresponding r-values 2, 
4 and 6 respectively. This example illustrates the big difficulty in getting these numbers by 
"brute force" for large values of N, (large means here iV > 10). 

2 Knuth uses the Christoffel symbols of first and second kind for Stirling numbers of first and second kind 
respectively. 
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Looking at (|^) we see that the number r is related by a particular partition of m, intercon- 
nected by another partition of n, both having the same height, i.e., the number r is exactly 
twice the height of these two partitions: 

















1 


1 


1 












1 


1 













1 

















1 


1 














1 


1 


1 



then, in this case r = 2h = 10. 

We define a "loop" associated to two partitions as the set of 2h strings obtained by 
matching these partitions, the first string regarded as following the last one, In the example 
we have: 



-00000 


— 111 — 00 




11 


HoHih-ooo 




11 




00 




111- 



For a fixed r, the set of all loops are obtained by matching the partitions of m with the 
partitions of n, both of height h = ^r. Fortunately this is done not by ordinary partitions 
but by partitions regarding the order, giving the numbers: 



M 



m 
h- 



1 



The total number of loops is given by: 



m 
h - 



1 



M 



n — 1 
h- 1 



Now we can get all sequences by making all N = m + n clockwise circular permutations on 
the digits in each loop, but in the process the h circular permutations displacing an even 
number of strings belong to different initial loops. Then, the total number of sequences is: 



Finally, taking 2h = 



N I'm — 1 
~h\h-l 
r our formula is 



—h 

nm 



t f m 



n 



T 



r = 2,4,6, 



~2 / \2 

The even index r runs from 2 to 2(mAn) and T 



1 1 
— + - 

m n 



(3) 



(4) 



mn 

T 



is symmetric in the upper indexes m, n. 
Table II. (T-numbers) 



h 


T 


N 


2 


3 


4 


5 


6 


7 


8 


9 


10 


m 


1 


1 


2 


1 


2 


1 


3 


2 


1 


3 


2 


1 


4 


3 


2 


1 


4 


3 


2 


1 


5 


4 


3 


2 


1 


1 


2 




2 


3 


4 


4 


5 


5 


6 


6 


6 


7 


7 


7 


8 


8 


8 


8 


9 


9 


9 


9 


10 


10 


10 


10 


10 


2 


4 






2 




5 




12 


9 




21 


14 




36 


32 


20 




54 


45 


27 




80 


75 


60 


35 




3 


6 














2 






7 






24 


16 






54 


30 






120 


100 


50 






4 


8 


























2 








9 








40 


25 








5 


10 










































2 
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From (|4]) we have the recurrence formula: 

T mn _ A m ~ k T )( n ~ k T ) -rmn 

and T™ n = number of clockwise permutations of the sequence 



000 0111 • • • 1 = m + n = N. 



In contrast with further generalizations, (j5|) shows the simplicity of these numbers: 

TT = N: 



2 

m — l)(n — I) 



TT = N , . o 

T mn = N (m-l)(7i-l) (m-2)(n-2) 
6 " 1-2 2-3 



(6) 



3. Properties of the Numbers T™ n 

At this point we can relate our problem to the probabilities. Since □ T™ n is the 
probability that a sequence have r jumps, normalisation gives: 

2(mAn) / 

y Tr = ( 

T=2 v 

i.e.: 

E" ,MM mn[N\ 



/i=i 



2 

E. / m\ m / 2m \ 



7l=l 



which is a well known formula (Knuth 0, p. 59])Q 



3 From 



Ji J \m — h / \ m 
(Knuth p. 58]), with m = r = s, we have also 



E([L*.J- r r. <ff) 



i x 7 

A method which we use later, consists of relating these formulas: 

i \ / o v 7 o v/ o v/ 1 



which gives (0"') from (f^). 
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For the ordinary moments of the above distribution we have: 
^VTT; 



V 



1 



m 



T' 



\m I z — ' mn \ h I \h I \mj mn \h I \h 



Obviously the probability of jumping isp = j^^jx, (probabilities of finding (00) and (11) are 
"1^-1) an( j ^fejjr respectively); then, according to binomial distribution, the probability 
of jumping exactly r times is: 



2pV(l-p)"- T (« (Z) TD (9) 



V r / \ m / 

(the factor 2 is due to the fact that jumps came by pairs), which shows how the theory of 
binary sequences would be complicated, even for physical problems, if one follows statistical 
formulations, because the exact formulas are much more simple than the approximate ones. 
As is shown in the following table, statistical approximation is good enough anyway: 

Table III. 



T 





2 


4 


6 


8 


10 







10 


80 


120 


40 


2 


Binomial 


0.15 


10.66 


77.71 


121.42 


40.65 


1.46 



Since both distributions are so close (at least in this example), we can expect a very good 
approximations identifying its moments. Call k r and (k) r the ordinary and factorial moments 
respectively. We have for a distribution po,Pi, ■ ■ ■ ,P2mA2n'- 

2mA2nAr 2mA2nAr 

F= J2 S rW)t= E S ' GW (1) (10) 

1=0 1=0 

the being the Stirling numbers of second kind and G^(l) the I derivative of the generating 
function G(t) at t = 1; (here we will be interested only in the case 2m A 2n A r = r). The 
generating function for Bernouilli's distribution is G(t) = (q + pt) N , so we have for ( jID[) : 

r 
£=0 

then, for r > 1, (for r = 1 we have the formula (0)): 

h){h)~{mAn) 2^N{N - 1) ^ S '- lA ™ (12) 

with 

A 1 =1 A' _f2mny- 1 (iV-2V 2 . 

mn ) mn I J ^jy- -|\£_2 > \ J" 
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As we will see in the sequel, approximation ([TJ) is extraordinary good for m — n, which 
gives: 



h=0 



2m N 
m I 2 T 



m " ( S i + mS 2 +m 2 S 3 2m ~ 2 
■i(2m-l) ^ b r-i+ mb r-i + m b -i 2 m-l 



P 2r.r-i 2m - 2 2m-r + 2 



2m — 1 2m — 1 

For discussion, we write explicitely ( |13"D from r = to r = 5: 

'2m N 
m 



(13) 



2m 
m 

2m 
m 

2m 
m 

2m 
m 

2m 
m 



— 1 m 



E fe0 

h=0 

-1 m 

E ftl 

/i=0 
m 

E* 2 

/i=0 
m 

E'* s 

h=0 

m 

E ,j4 

h=0 
m 

E'* 5 



h=0 

-1 m 



h=0 
-1 m 



h=0 
-1 m 



h=0 
-1 m 



1 



(0,0) 



h=0 



m, 
~2 



m 



2{2m - 1) 
m 

4(2m - 1) 
m 

8(2m - 1) 
m 

16(2m- 1) 



[l+m) 



1 + 3m + m" 



,2m -2 



1 + 7m + Qm z 



2m - 1 
2m - 2 



©2) 
©3) 
©4) 



2m - 1 



+ m' 



(2m-2)(2m-3) 
(2m - l) 2 



(0,5) 

We know that formulas (0 0) and QE| 1) are exact (<g"), <g)) and we will show that (0, 2) 
and fll3| , 3) are exact too and also why identity fails after r = 3. 

Proof. 

2 m / / \ \ 2 



E'» 2 7 =E * 



m / i \ 2 

m — 1 



m 



E 



m 



! \ \ / / ! 

2 m(2m — m) f 2m N 



/i-l 



m 



2m -2 
m — 1 



2m(2m-l)\m/ 2(2m - 1) \ 



m 3 / 2m 



m 



In the general case: 



m / \ 2 m / \ 2 r 

E-I =E(-*)'Q =E< 

i v 7 o v 7 fc=0 



m 



•-fc 



E*" 

>/i=0 



then, we have h r (™Y by recurrence onfo/ /or og?g? values of r, because when k = r = 2 the 
last term cancels with the left hand side, giving only an identity between moments of lower 
orders. Only the case r = 2 has escaped from that inconvenient (due to the square in (™)) 
and provides for ( [IB] , 3). Thus, for r = 3, we have: 
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m / __\ 2 771 / \ 2 N 

3 I 



rn / \ z 771 / \ z 7n / \ z 771 

E ft) -*»'E*ft) ^E» 2 ft -E* 



i.e.: 

'm\ 1 / o 9 m m 3 \ m 3 (m+l) 

' - 'm 3 -3m 2 - + 3m ' - v 



h) 2\ 2 2(2m-l)/ 4(2m-l)' 

Recurrence fails for /i 4 (™) 2 and so for higher orders. □ 

Formula fll3[ 4) gives the following approximations. For m = 2: 30/9 ~ 29/9; for r = 3: 
11.70 « 11.61; for m = 10: 827.40 « 827.22; and for m = 15; 3829.74 « 3829.48. Formula 
(0,5) gives, for m = 10: 4895.51 « 4891.65. Then, the distributions ( 2 ™) _1 T™ n and 
2( 2 ™)p T (l — p) N ~ T have the same variance; higher order moments are almost equal in both 
distributions. 

Let us look for these relations in case m ^ n: Supose m < n (the case m > n is the same); 
according to (^") we have: 



/i=0 



h I \h I ^— ' \ h J \n — h J V n 



then: 

TO. 

e^i:)(:)=e(m:))(m:))=™e 



mAn / \/ \ 77*//\\//\\ 77i— 1 / .. , , , 

m\ / n\ s-^ / . / m\\ I n\\ im — 1 \ /n — 1 



"I— 1 / 1\ / 1\ / n\ 2 2 /AT 

m— l\/n — 1\ I m + n — 2\ mn f JM 



E/ 7/4 — 1 \ / 74 — 1\ / 7/4 -|- ft — 

I Zi A /» =m H m-1 



o 

which is the formula ([12]) for r = 2. For r = 3, formula (|P2| ) gives: 



\ m-1 / iV(iV - 1) \m 



E n 3 fm\ fn\ ^ f N \ / 1 m 2 n 2 m 3 n 3 

/i=l 



/i/VV \mAnj \2iV(iV-l) iV 2 (iV-l) 



which is only approximate; for m = 4, n = 2, (N = 6) gives 56 ~ 58.67. Of course the 
approximation will be much better for higher N. Then, for m ^ n, only the mean values 
are equal in the above distributions, though higher order moments are very close in both 
distributions. 

Let us consider the set of all 2^ sequences of length N. One may expect here complete 
agreement with probabilities, i.e. jumps must follow the binomial distribution with p — |. 
The next proof shows that this is the case: 

Proof. From 

y-v ( a — c \ (h ~i r\ ( a : h ; ] 

0<c<a 



3 / V & / \J + k + 1 
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(Knuth |], p. 58]), we have (b = 0, j = k): 



El a — c\ I c\ f a + 1 
V k ) \k) = \2k + 1 

0<c<a v / \ / \ 

Then the exact number of sequences with r jumps ((01) or (10)) in all 2 N sequences is, 
according formula (H): 

T iv-n,n N^((N-2)-(n- l)\ (n - V 



n=0 



^ T " h ^ \ h-i J U-i 



h\2h-lj \2h \t 



(14) 



So, the probability of jumping r times: ^ n T™"/2 7V = (^)/2 N 1 equals that given by 
Bernoulli's distribution <g) with p = \: 2(f) (|) r (%) N ~ T , (r = 0, 2, . . . , JV (or iV-1)). □ 

Application to the Ising model of ferromagnetism. Let us consider a one-dimensional 
ferro magnet in absence of external field and put v = J / kT (J = interaccion energy of antipar- 
allel spins (—J for parallel); k is the Boltzmann constant and T the absolute temperature). 
The energy of a given distribution of iV spins is J(—(N — r) + r) = J(—N + 2r), because 
there are r couplings between antiparallel spins. The partition function for a fixed number 
oi N — n and n spins up and spins down respectively, is: 

% = ^ ■jN-n,n e -{-N+2T)v 

T 

(this is also the partition function corresponding to binary alloys (ions A and B instead 
of spins up and down)). Fortunately (for the Physicists) one must take the sum of all 
configurations which gives: 



Z = Y J Y. T?- n > n e-(~ N+2T » = 2e Nv ^ ( N \ (e^T 



= e Nv {{\ + e~ 2u ) N + (1 - e~ 2u ) N ) = (2 cosh v) N + (2 sinh u) N . 

One wonders why some authors give the solution which contains only the hyperbolic cosinus 
(see, for example, [0, ch. 5]; sinh appears in more elaborate discussion || ch. 4, p. 121]) 
though both solutions are equal for iV — ► oo. We think their countings are not accurate 
enough, for they take the following recurrence: Z N = (2 cosh u)Z N _ 1 on grounds that a 
new spin provides an ammount of ±J in the energy, i.e. forgetting a border effect. In fact, 
consider three spins: A up, B down and C up and add a new spin D; the old energy is 
2 J — J = J, and the new energy is for D up or 4J for D down and not J ± 2 J given by 
the above rule. 

4. Asymptotic Expansions 

After iV = 8 we can use the asymptotic approximation for T™ n , with the help of Stirling's 
formula which gives: 



\ 2 m+l (15) 



h ) v / 2vrm 
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(|D)p is surprisely accurate even for low values of h (< m); m — 10, h — 3 gives: 120 = ( 3 °) ~ 
116.1. (For h = \m we have a/^v = ~/^{in) 1 which is, of course, the Wallis' formula 
because (i^J = 4(/^)^| gives n N+2 = t^nj^j^- and since 7r 2 = 2 = 4, we have: 



7T4 



2-2 4 
TT3'3 ;7r6 



2 ■ 



2-2 4-4 



6 

1^3 ' 3~5 ' 5 ; "- ;7T ° 



2-2 4-4 6-6 



7T 



1-3 3-5 5-7 

We can expect a good approximation for T™ n when A" > 10; in fact the following formula 
gives, for iV = 10 a result very close to the Table [|: 



1 h ~ 



i i „-2?t 2 fi + l)+4ft.-i(m+n) 



h - + - 



v m n j 



or: 



-rmn 
1 /i 



2Ne 



ir(mn) 3 / 2 



m n j 2-RJmn 



he~ 4 ^r h - 



-- + log2 « 0.193147). 



In terms of r's 



(16) 



(17) 



For A" = 10, m = n = 5: T^ 5 = 0.3514 /iexp ( — 4(^- — /i)) and we have the following results 
as is shown in Figure [I]. 

Figure 1. Asymptotic expansion for T^ 5 
Tf = 0.3514 xexp (-4(f - x)) 



150 



100 



50 - 
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Table IV. 



T 2 


4 


6 


8 


10 


12 


Tf 10 


80 


120 


40 


2 





Asymptotic 8.62 


80.40 


127.95 


34.44 


1.76 


0.02 



At the origin, the derivative of (|i~6|) is not zero but 2(7rny/ran)~ 1 e aN and the curve takes 
though very small, negative values for x < 0. 

On considering the set of all 2 N sequences, the distribution of jumps, for a large N is 



with (|rj) and (||): 

Here the gaussian has nothing to do with statistical distributions, because in this context 
(|P8]) gives fluctuations, respect to the exact number of jumps, due only to the fact that N is 



not large enough. Of course ( |i8|) becomes a normal probability distribution with standard 



error ±v iV when the sequences of N digits are taken from a table of binary random digits. 

5. The Numbers T^"(00) and T^ n (ll). 
Generalized Fibonacci Numbers 

Let U be an ordered sequence of r zeros and s ones, \U\ = r + s < N; being r the 
amount of (01)'s or (10) 's in the sequence Si, our general problem will consist in finding the 
number of times, h = h(U), a fixed string U is contained in Si in the following sense: If 
Si = a.\ (%2 «3 • ■ • «at, {a.j = 0, 1), take the N strings: 

(aii a 2 a 3 ■ ■ ■ a r+s ), (a 2 a 3 a 4 ■ • • a r+s+ i), . . . , 

(a!jv_i a N a x ■■■ a r+s _ 2 ), («at ai «2 • • • «r+ s -i), 

then h is the number of these strings having the form U. For instance, take in (0) U = (1100), 
then h = 3. We call T™ n U the number of sequences which have h times the string U. 

Now, since after any (01) in a sequence S sooner or later another appears (Ising model!), 
then h(10) = h(01) = |r, i.e., 

T^ n (01) = Tr(10) = T™ n . (19) 

The numbers T™ n (00) and T™ n (ll) are found by a simple consideration: From all iV pairs 
succesively taken in S, m = the number of zero digits = the number h + |r, of them having 
the form (00) or (01), so the r corresponding to S is 2(m — h), and we have: 

TT(00) = T^_ h) , similarly: TT(H) = T^_ h) (20) 

then: 

Tr( oo) = n'^Ji h)( n X miD = { m ,) (") (2 i) 

ran \h J \m - hj 11 ran \m - hj \h J 

which gives the obvious identity: T^ n (00) = T^ n (ll). 
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Normalization. Since the only values of h for the T™ n (00) are between m and m — (mAn), 
we have, (the same for T™ n (ll)): 

mAn mAn mAn mAn 

£ T r(oo) = £ Ts^ h) = x: "Hy" =£r = 

m—h=0 m—h=0 j=0 5 r= 

Let F*(iV; r, /i) be the number of subsets of X = {1, 2, 3, ... , N} that do contain exactly 
h times r consecutive integers modulo N; as in Berge @, p. 31] let us associate with each 
subset Y C X a "word" S = a±, 0:2, 03, . . . , Qn, where a, = if i Y and ctj = 1 if i £ Y. 
The mapping between subsets Y and words S being bijective, instead of counting subsets, 
we can count the number of the sequences S in which the string (11 • • • 1) is repeated exactly 
h times. We call Generalized corrected Fibonacci numbers the number of these subsets, i.e.: 

F*(N;r,h) = J2<~ n ' n (^- : -V (22) 

n 

then: 

n 7i=/l+l V y V / \ / 

Example. X = {1, 2, 3, 4}, r = 2, h = 2. In this case, in ( [Z3"D the only index is n = 3, then 
F*(4;2,2) = 4. The subsets are: {1, 2, 3}, {2, 3, 4}, {3, 4, 1}, {2, 4, 1}. For h = one obtains 
the corrected Fibonacci numbers: (see Berge, @, p. 32]) 

F* (J V;2,0) = F* = £— ( j. 

71=1 V 7 

Next, to each of Q = (*) sequences we will associate a fa/pe, which is the pair of partitions 

with height h, i.e.: X^i a « = Si a i = ^5 corresponding to the tableaux of 0's and l's, 
respectively, disregarding the order of the blocks, i.e.: in descending order of the length of 
blocks, as is customary. Two types are equal iff the have both Young's tableaux equal. For 
instance, for m = 4, n = 3, the 35 sequences of the Table | have the following four types: 

1 = {(4),(3)}, 2 = {(3,1), (2,1)}, 3 = {(2 2 ),(2,1)}, 4 = {(2, l 2 ), (l 3 )}. 

Thus, 

type{0011100} = 1, type{0110100} = type{0001011} = 2, etc. 
In total, for a fixed m, n, there are the number 

mAn 

£ Pk m P n ( 24 ) 
h=l 

of types. The great number of sequences compared to the small number of types is due to 
the M-partitions; in fact, the number of sequences per type is: 

Nh\(h-1)\ 
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because the M-partitions are found making the sum: h\ ^2 ai \ a l\---a \ °^ a ^ partitions with 
same h. This clarifies the meaning of the formulas: 

TT(01) = TtlO) = ^M h m M h n , 
N 

T™ n (00) = — M™ M™ , (26) 
n 

tt(h) = -KT h W h - 

n 

6. The Numbers T™ n U for \U\ = 3 

The number of formulas can be reduced by some general considerations; for example, 
if U' is the string obtained by interchanging in U the 0's and the l's we have: T™ n U' = 

k 

T^U. Furthermore, since sequences are unbounded, after any string (011 ■ ■ • 1) in the 

k k 

sequence sooner or later another zero appears, then T™ n (ll • • • 10) = T™ n (011 • • • 1); also, 

TX"(00 ••• 01) = T™ n (10 ••• 00), so Tr(H ■ 10) = T" m (10 ■ •• 00), etc. For U = 3 we 
have the following identities 

T™"(001) = T™ n (i00) = T^llO) = T^ tm (011); 

TT(lOl) = TT(OIO); 
T™ n (lll) = T" m (000). 
Thus we only need formulas for the three expressions: 

T™ n (001), "nT(lOl), T^(OOO). (27) 

Let P^be the number of partitions, disregarding order, of weight m and dimension h; as 
is well known, the process of deleting the first column on all associated Young's Tableaux 
gives: 

s 

P h m = J2 p£ m-k {s = h A (m — h)). (28) 
Let us try a similar expansion for M-partitions: 

s 

< = j>Mj>_ fc , (29) 



i=i 



formula (0) gives: 



t(';)M,u^t<';)('7 ft )^C:)^f, (so. 

Then the coefficients a# in the formal expansion ( p5|) are equal to Q). We shall prove that 
this expansion is not only formal which proves the following lemma: Define c* fc = ( t ~ k :> ) 
for i ^ j and c\ k = iSoj., so that (0) has the form: c* fc = M]. We have, 

Lemma I^ 1 ^. The number of M-partitions of m in h parts is but after deleting the 
first column in all Young-De Moivre Tableaux we have M-partitions of dimensions i = 
0, 1, 2, . . . , h A (m — h), which number is c^. 
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In order that the counting makes sense we must define c™ = 1. Further, it is important to 
take in consideration the fact that deleting the first column is made for counting purposes 
only. Consider, for instance, this process for m = 6, h = 2; we have the (J following 
tableaux: 



I 



m 



-20 







m 



21 



r 6 

c 22 



then c®i = 2 means that the tableau I I I I I shall be counted twice because it has different 



sites with respect to the first column. 

Proof of Lemma ^ . Before deleting the first column, the number of tableaux which will 
give tableaux of dimension I after the cut, is (Figure 0): 



Figure 2. 



7 



h 



«3 



a 2 



h-e 



r \- h\ h\ I 



(h - £)\ a 2 \ a 3 \ ■ ■ ■ (h-£)\ ^a 2 \a 3 \--- 
where the sum is extended over all tableaux. After deleting the first column, these tableaux 
will give M e m _ h tableaux of dimension £, which number is: 



then: 



. . m—h) 

a 2 \ a 3 \ 



(h-t)\ £\ ~ 1 " ' m ~ h ~ ht 
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Due to importance of coefficients c* fc we give the tables (for matrices of low k) in Appendix 
A. It is immediate to find a recurrence relation: 

a U ~ - k) 



c j,k+l 



kJkTT) • (31) 

Now it is a simple task to find T™ n (001); as in the case of T™"(01) we match tableaux 
of zeros with tableaux of ones of same height h. Since strings (001) are originated from the 
rows of 0's tableaux which remain non empty after deleting the first column, here in the 
first formula fl26|) instead of we put c^, which gives a number t of (001) strings; these 
t strings are a subset of the h strings (01), {h > £), because both (001) and (101) give (01). 
Then we have an intermediate formula for (001): 

. The number of sequences in which the string (01) appears h times and the string (001) 
appears £ times is given by: 

>(01;001) = ^CM£ (32) 



(Note that: J2e T S"(01; 001) = f M£ J2 e which is formula @). Thus T^(OOl) is found 
by disregarding h in (|32]), i.e.: 



T™ n (001) 



N 



n 



"EfM5= T E<Og- (^) 

h h V 7 

Let us find the exact limits on the sum; obviously we have: £ = 0,1,2,..., [y] A n, and 
the inferior limit corresponds to h = i. Condition h < m,n must be restricted because in 

to k = i — j, then h < m — £,n and finally: 



c l j k the low index k runs from k 



-rmn 
1 I 



(001) 



N 
n 



c 



(m-£)An 

E 

h=e 

(m—£)An 



III! 



(33b) 



N ^ 
n ^ 



m 



h 



h\ ( m — h 



In that follows some binomial coefficients will be taken outside the sum by help of the 
formula: Q Q = g) { a J) , then « is: 

TT(001) = -( ] Yl )[„_,)> (n) =1 - ( 33c ) 



Let us take an example which can be checked by hand on the 70 sequences of 4 zeros and 
4 ones. From the tables of Appendix A we have: 







h 





1 


2 


3 


4 















1 




1 








2 




2 


1 






3 




3 








4 


4 











then T 4 ^(01;001) = -c; 







h 





1 


2 


3 


4 















1 




8 








2 




24 


12 






3 




24 








4 


2 
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The sum rows and columns gives Tf (01) and Tf (001) respectively: 





i 




h 





1 


2 


3 


4 


TT(oi) 


















1 




8 








8 


2 




24 


12 






36 


3 




24 








24 


4 


2 










2 


17(001) 


2 


56 


12 


9 








Thus: Tf (001) = 2, Tf (001) = 56, Tf (001) = 12. The first corresponds to two sequences: 
{01010101}, {10101010}. 

Taking in ( |33b| ) the sum over m, we have also the distribution of the string (001) for all 
2 N binary sequences of cardinality iV: 

Tfcoi,-* ± T {N j ){m _ h) Q( m - h )( N r) 

m =e+i h=l v /v J \ / \ / \ / 

~J^( 1-1 £-1 )\ h-£ )' 

Turning to the second expression in (p7|) , we can see in Figure |] that to each tableau 
of dimension I obtained by deleting the first column, there corresponds the bottom h — I 
positions of h which are responsible for the strings (101), then: 

. The number of sequences in which the string (01) appears h times and the string (101) 
appears £ times is given by: 

N i N /n\ 

T£ B (01; 101) = T™_,(01; 101) = ^ h _ t) M h n = -c% h _ t) . (35) 

Moreover: 

T(ioi) = ^E c M^)(l)- ( 36a ) 

Here, the index £ runs from to m A (n + 5 mn — 1) because the extreme cases have the form 
(for m > n, m < n and m = n): 

{000- • -0101 •• -1010- • -000}, {111 •••10101 •••101- ••111} and {10101 ••• 1010}. 

In the sum, index h goes from h = £ to h = [^^j (because h = £ = m gives c™ = 1 and 
m — h = h — £ for the high index), then: 
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The number of sequences of m 0's and n l's which have i strings (101) is: 



h=l 

[2^]Af» 

N (n\ sr-^ f m — h — 1\ (n 



(36b) 



There remains only the last expression in ([27]) ; however we must first see how the numbers 
T(000) are related to string (0001) and so we must pospone its derivation until after the case 
\U\ =4. To justify this delay and find a way to obtain T(000) we will see that, though more 
complicated than (|20|), formula (|21"| ) for T(00) can be found also in terms of strings in which 
\U\ =3. Let g be a fixed number of strings (00), it is easy to see (Figure |2|) that m — h = g; 
i.e., there are as many (00) as zeros in the subtableau of length £. After deleting the first 
column of length h = m — g there appears a number c™_ g £ of such subtableaux, in total an 
amount of gc r £ l _ g l zeros, then: 

m ~9 AT m ~3 AT D / \ / \ / 

N m K>m-Q \- N * ( m ~ g\ i g\ m - 9 1 n 



~[ m — 9 m ~ 9 9 \ " J \ 

= %>n- g) l?( m ~i g ) © = ^ <m " 9> (r) 



n \m — g 



£=0 

The last derivation is due to: 



m \ I n 
m — g 



m-g / \ / \ 9 / \ / \ (m-s)Ag 

?<9(v)-?<9(v)- ?'" wm 



allowing to use formula (|7|). 



7. The Numbers Tf n (00 • - oi). 
Generalized Kaplansky Lemma 



Define: 



jA(i-j-k) 



A A-i-Z 



Jfc " ^ i-j-f\fj\kj\ k 



(37a) 



and: 



then, 
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Lemma I^ 2 ^. After deleting the first and second columns in all Young-De Moivre tableaux of 
weight m and dimension h one obtains M-partitions of dimensions £ = 0, 1, 2, ... , h/\{m—2h) 
which number is c'^. 



The proof is like 1^, (see Figure 



Figure 3. 




The desired number is, (the sum without indexes is extended of all M-partitions): 

Af= ^^(h-£')\(£'-£)\ ' a 3 \aA ■■■ = ^ ^QQ^ a 3 !a 4 ! •■■ 
on the other hand: 



E 



n 



Then: 



h/\(m-h~e) 



«E 



1 



03! CK4! ■ 



M 



m—h—l 1 " 



N 



E 

i'=l 



M 



m.-h-e 



M 1 



(^0); 



the upper limit is due to the number of elements m — h — I' in the subtableau which must 
be at least equal to t. In the case £ = the symbol M° m _ h __ e , is meaningless. In this case 
however we observe that the condition m — h — £' = gives M = c'™ = ( \) • 
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. The number of sequences with m zeros and n ones in which the string (01) appears h times, 
string (001) appears £' times and (0001) £ times is: 



x N fh\ fh- A fm-h- £' - 1\ fn 
N fn\ fn -l\fn-l'\fm-h-t'-l 



(38a) 



And in case m < 2h: 



T^,o(01; 001; 0001) = f » J (" " ™ ; A ) , (38b) 

which, of course, coincides with T^,_„(01; 001) = fcj m _ ft g), because cj^ = (,^), 
Though binomial coefficients make elementary computations, it is always convenient to give 
fundamental formulas in terms of c's because its matrices can be tabulated independently of 
particular cases. We have for T(0001)'s: 

Tr(oroi) = ^xy«M* 

{m-2t)An , „ , w . (39a) 



TT(0001) = - £ ( 

Hf] 



?C) s ?(::;:)(•;:;)("-;:■'-') «« 

(39b) 



nhira—n) 



n \ I n — m + h 
m — hj \ 2h — m 



For the general case we define: 



^ ( fj N /' /'■ 



(40a) 



i-j-fi-f2 f, 

k 

k ) j , /; , ; ,, V/i - iv U - / 3 / " U-i - /J U - 
< ~ i ~ /l ~/-~i""~ / "~ 1 ) (M0) ' 

,.< "■)' - J ^3>h>h>->fs (/1-/2) (/2-/ 3 ) ' ' ' (/.-1-/J (/J ' " ( " ' ' ) - / - ,101,1 

Cj0 ~ i if<>(s + i)j; 1 j 



then, 



Lemma I. deleting the s + 1 columns in all Young-De Moivre tableaux of weight m 

and dimension h one obtain M-partitions of dimensions £ — 0, 1, 2, . . . , h A (m — (s + 
which number is c^ m . (This generalizes I* 1 ) and 1^.) 
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So, we have: 



T^»...^)(01; 001; • ■ ■ ; 00 ( - • - } 01) = jc^ m M h n , 



c 



" ' un being the c^j m without summation on £', £", . . . , £^ s \ And: 



// \ \ Oil Dm 



h=t e'>£">--->e(")>h 



n 



'm - h - £' - £" £ (s) - f 

h-£j\ £-1 

(41) 

Formula (|4lD generalizes Kaplansky's lemma (for a redefinition of Kaplansky's lemma see 
Comptet p. 35]). In particular, taking h = £' = £" = ■■■ = $> 8 > = n, we have: 

/AT . N fm — n — sn — l\ N [N — pn\ 

g n (N,p) = -( )= , Kaplansky 

n \ n — 1 / N — pn\ n J 

being p = s + 2 the least length of blocks which separate the ones. 

8. The Numbers T™"(00 - - 00). 
Fibonacci F*(N;r,h) Numbers 

Our final problem will consist of looking at the weights g of the subtableaux, instead of 
the dimension £. We will call c^ mg the number of subtableaux of weight g, which remain 
after deleting s + 1 succesive columns, starting from the first, in all the Young-De Moivre 
tableaux of weight m (> g) and dimension h. 

For s = 0, in Figure ^| we can see that the weight of the remaining subtableaux is m — h, 
then (II* 1 )): 

r ™ a _ST r m _ { ( m ~ 9\ f9\ _ m ~ 9 fm\ 

h ~h 9 ' l ~o~a\ * )v)~ - \g) 



m — g ( m \ fm — 1 



(42) 



c mg 



m \m — g J \ g 

so c™ 9 = c m9 does not depend on h, i.e.: Counting shall be made on all tableaux of weight 
m. Though indirectly, it is thanks to numbers c^ g that T(00) was found and we may now 
write formula (|l]) as ^c mg M™- g . 

For s — 1, it is easy to see that the problem of counting strings (000) is given also by the 
formula ( p9|) but collecting the terms with the same g, i.e., from (|37| ) we obtain c'J taking 
in the sum only the term / = i — j — g: 
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We have (see Figure 



img 



E E 

i m—h—t'=g 
hAgA(m—h—g) 



hAg 



M 



-h-e - 



E 



t 



£=0 



h 



in — g — h — I 



h-e 

m — g — h — £ 
(g + 0); 



(43a) 



JmO 
C h 



hO 



L-h) '*m<2h, 
if m > 2h. 



(43b) 



Then we have: 



Lemma II^ 2 ). After deleting the first and second column in all Young-De Moivre tableaux 
of weight m and dimension h one obtains M-partitions of weights g = V (m — 2h), 1 V (m — 
2h), . . . , m — h — 1 which number is c 1 ™ 9 . 

In order to clarify ambiguity, let us make the process for m = 7, h = 3: 



J 73 



n 



5 



m 



J 71 



m 

; — 9 (two types 

m 



i 



J72 



In c'g 2 tableaux, first and fourth give the same position for the subtableau of weight 2, but 
since the second column is different that will give different matchings; for example {n = 9), 



































1 


1 


1 


1 1 















1 


1 


1 


1 1 


1 


1 


1 




and 









1 


1 


1 




1 


















1 









give different sequences, 

{0000111110011101} and {0000111110111001}. 

Now, at last, the coefficients c'™ 9 provide for the numbers T(000) (and so for T(lll)). In 
fact, there are as many strings (000) as zeros (g) in the subtableaux of weight g, so that: 

IT (000) = £ £ cT = £ £ c'r (f) . 

h h \ ' 

Since the matrices c'™ 9 (g fixed) can be tabulated very easily, as in Appendix A, it is 
preferable to work with them; nonetheless if we write explicitly all binomial coefficients, we 
have: 



(44) 



-rmn 

9 



(000) 



N 



n 



EE 



?? 



n-£ 
m — h — g 



n — m + h + g 
2h + g — m 



(g ± 0) (45a) 



22 Structure of Binary Sequences © J. Tharrats 1982-2001 Printed 2001-0924 math. CO/0109182 



h h / ' 

Let us find a complete example with m = 5, n = 3, (N — 8) which can be easily checked 
by hand on the 56 sequences. Appendix A provide for all coefficients: 

T?(000) = 5ef (3) = IsQ) = 8; 

Tf (000) = § (cf g) + cf g) ) = § (2 g) + sg) ) = 24; 
Tf(000) = ^cf Q=16; 
Tf(000) = ^fQ=8; 

Tf (111) = Tf (000) = ? (c'fQ + c'fQ) = 48; 

Tf (111) = |c'f = 8; 
Tf (001) = 0; 

Tf (001) = I (cf Q + 4 Q + cf Q ) = 32; etc. 

Tf(011) = Tf(100)=Tf(001) = ^c? Q+4Q+4Q) =16; etc. 
The following table gives the full results: 



Table V. 



T 


(111) 


(110) 


(101) 


(100) 


(011) 


(010) 


(001) 


(000) 


T 53 
1 


48 


16 


24 





16 


8 





8 


-1-53 
1 1 


8 


40 


24 


32 


40 


32 


32 


24 


-1-53 
1 2 








8 


24 








24 


16 


-1-53 
1 3 

















16 





8 



In general, for s > 1, take i — j — f\ — f 2 — ■ ■ ■ — f s = g and define: 



and: 



"'"--f E k (i)(^)-( f ^(l)- U>f<>k-,9*0), (46a) 



9 k=l h>h>->fs V./i/ VAv \J> J v< 

/i+/2H \-f s =i-j-g 



w, °- E (j) ©-ftO' « (46b) 



h>h>->fs 
h+h+-+fs=i-j 
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(From hereon in order to calculate the c's one has to deal with the number of partitions 
P|_j, P|_,-_ ff , etc. which complicates the formulas). Then we have: 

Lemma II. After deleting the first s + 1 columns in all Young-De Moivre tableaux of weight 
m and dimension h one obtains M-partitions of weights g = V (m — (s + l)h), lV(m- 
(s + l)h), . . . ,m — h — s which number is c^ mg . (This generalizes IlM and 11^.) 

This lemma provides for U = (00 ■ • • 0) or U — (11 • • • 1): 

T r( 00 't 2 0) = Y, = £ £ c<«> m » (") , (47a) 

h h \ ' 

and: 

Tnil- +2 l) = ^£<Wf)- (47b) 
h ^ ' 

So, the generalized Fibonacci numbers are: 

(s~2)ng , N \ 

f-(jw) = "EE^(V > < 48 > 

n h x 7 

For example, the number of subsets of X = {1, 2, 3, ... , A^} that do not contain three con- 
secutive integers modulo iV is: 



N N-n , nQ , Ar s n +W\< N N-n I h \ (N-n\ 

^N-.m^N^j:^ r) =N S E%^ < 49 > 

Example. In ( j49|) take A^ = 5, we have: 



F*(5;3,0) = 5^0 + + I( 3 ^)n + i 



= 20. 

The twenty subsets are: 

{1, 2, 4}, {1, 3, 4}, {1, 3, 5}, {2, 3, 5}, {2, 4, 5}, 

{1, 2}, {1, 3}, {1, 4}, {1, 5}, {2, 3}, {2, 4}, {2, 5}, {3, 4}, {3, 5}, 

{1},{2},{3},{4},{5}. 

Appendix A 

The recurrence relation (|31~D gives the matrices c* fc , we have, {i rows, j columns, /c fixed): 
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10 



11 



12 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


























2 


1 
























3 


1 


2 






















4 


1 


2 


3 




















5 




2 


3 


4 


















6 




2 


3 


4 


5 
















7 




2 


3 


4 


5 


6 














8 




2 


3 


4 


5 


6 


7 












9 




2 


3 


4 


5 


6 


7 


8 










10 




2 


3 


4 


5 


6 


7 


8 


9 








11 




2 


3 


4 


5 


6 


7 


8 


9 


10 






12 




2 


3 


4 


5 


6 


7 


8 


9 


10 


11 






1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


























2 


























3 
4 




1 






















5 




2 


3 




















6 




3 


6 


6 


















7 




4 


9 


12 


10 
















8 




5 


12 


18 


20 


15 














9 




6 


15 


24 


30 


30 


21 












10 




7 


18 


30 


40 


45 


42 


28 










11 




8 


21 


36 


50 


60 


63 


56 


36 








12 




9 


24 


42 


60 


75 


84 


84 


72 


45 
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10 11 12 



C }3 = 



1 
2 
3 
4 
5 
G 
7 
8 
9 

10 
11 
12 



1 












3 


4 










6 


12 


10 








10 


24 


30 


20 






15 


40 


60 


60 


35 




21 


60 


100 


120 


105 


56 


28 


84 


150 


200 


210 


168 


3 


4 


5 


6 


7 


8 



84 



10 



11 



12 



C j4 - 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 



1 

4 5 

10 20 15 

20 50 60 35 

35 100 150 140 



70 



9 10 11 



12 



u j5 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 



1 

5 6 
15 30 



21 



And the finite matrices c l jk (i fixed, j rows, k columns) are: 
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jk 








1 


2 


•> 

3 


4 


r 











1 


2 


3 


4 


5 


1 




1 












1 




1 










2 




2 












2 




2 


1 








3 


3 














3 




3 










4 
















4 


4 












R 
u 














c 4 , - 


R 

K 
u 














7 

a 
















7 
« 














O 

9 
















g 














10 
















10 
















U 


1 


o 
z 


Q 

o 


4 


r 








U 


i 


o 
z 


Q 

o 


4 


r 




1 




1 












1 




1 










2 




2 


2 










2 




2 


3 








3 




3 


3 










3 




3 


5 


1 






4 




4 












4 




4 


6 








5 


5 












c 6 , = 


5 




5 










R 
u 














a 

\j 


u 












7 

a 
















7 

a 














O 

9 
















O 

9 














10 
















10 
















U 


1 


I 


Q 

o 


4 









u 


1 




Q 

o 


4 





1 




1 












1 




1 










2 




2 


4 










2 




2 


5 








3 




3 


9 


3 








3 




3 


12 


6 






4 




4 


12 


4 








4 




4 


18 


12 


1 




R 




R 


1 f) 








c 8 , - 


R 




R 


20 


i n 






u 




fi 
u 










a 




a 

V) 


1 r i 








7 

a 


7 














7 

a 


Q 


7 










O 

9 
















O 

9 


o 












10 
















10 



















1 


2 


3 


4 


5 









1 


2 


3 


4 


5 


1 




1 












1 




1 










2 




2 


G 










2 




2 


7 








3 




3 


15 


10 








3 




3 


18 


15 






4 




4 


24 


24 


4 




„10 _ 
c jk - 


4 




4 


30 


40 


10 




5 




5 


30 


30 


5 




5 




■5 


40 


GO 


20 


1 


6 




6 


30 


20 






6 




6 


45 


GO 


15 




7 

8 




7 

8 


21 










7 

8 




7 

8 


42 
28 


35 






9 


9 














9 




9 










10 
















10 


10 













Let us see how the matrices c'* fc can be obtained from the matrices c* fe . Formula (|37|) can 



be written: 



7* 



i.e., C' k = C k V 



(1,A) 
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where (C' k )ij = c' 1 ^ 3 , (C ) = c , and (P) rs = Q is the Pascal matrix, so (r rows, s columns): 



(V) ra = 





1 


2 


3 


4 


5 


6 


7 ••• 


1 


1 


2 


3 


4 


5 


6 


7 ••• 


2 




1 


3 


6 


10 


15 


21 ••• 


3 






1 


4 


10 


20 


35 ••• 


4 








1 


5 


16 


35 ••• 


5 










1 


6 


21 ••• 



We have succesively (i rows, j columns): 

12 3 4 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 



1 



9 10 11 



10 



11 





1 


2 


3 


4 


5 


6 


7 


8 


9 


1 




















2 




















3 


1 


















4 


1 


2 
















5 


1 


4 


3 














6 


1 


4 


9 


4 












7 


1 


4 


12 


16 


5 










8 


1 


4 


12 


28 


25 


6 








9 


1 


4 


12 


32 


55 


36 


7 






10 


1 


4 


12 


32 


75 


96 


49 


8 




11 


1 


4 


12 


32 


80 


156 


154 


64 


9 


12 


1 


4 


12 


32 


80 


186 


294 


232 


81 



10 11 



10 
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9 10 11 



J i 
C j2 



-j3 



Finally, formula 



1 

2 
3 
4 
5 
G 
7 
8 
9 
10 
11 
12 



1 












2 


3 










3 


9 


G 








4 


15 


24 


10 






5 


21 


48 


50 


15 




G 


27 


72 


120 


90 


21 


7 


33 


9G 


200 


255 


147 


2 


3 


4 


5 


6 


7 



28 



9 10 11 



1 






2 






3 






4 






5 






6 






7 






8 
9 


1 




10 


3 


4 


11 


6 


1G 


12 


10 


3G 



■kj 



10 
50 

can be written 

J-f 
j-(k + 2f)-(j 



E 



20 



*0 



which gives successively (i rows, j columns): 



J 

3-f 



k-r 
f-ij' 



0j 





1 


2 


3 


4 


5 


6 


7 


8 


9 




1 


2 


3 


4 


5 


6 


1 


1 


















1 














2 


1 


1 
















2 














3 




2 


1 














3 


1 












4 




1 


3 


1 












4 




2 










5 






3 


4 


1 










5 




2 


3 








6 






1 


6 


5 


1 








6 






G 


4 






7 
8 








4 
1 


10 
10 


6 
15 


1 
7 


1 




7 
8 






3 


12 
12 


5 
20 


6 


9 










5 


20 


21 


8 


1 


9 








4 


30 


30 


10 










1 


15 


35 


28 


9 


10 










20 


GO 



8 9 
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-2j 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



4 

18 5 

30 50 6 

10 60 45 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
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2 
4 



3 

12 4 
20 24 



40 40 



This work was partly done while in New York University on sabbatical leave. I thank Prof. 
Martin Pope of the Radiation and Solid State Laboratory and Prof. Peter Lax of the Courant 
Institute for their hospitality, and Prof. Frank Anger of the Department of Mathematics of 
the University of Puerto Rico at Rio Piedras, for useful comments on the manuscript. 

References 

[1] Andre, D. Notice sur les travaux scientifiques. Gautier Villars, Paris, 1910. 

[2] Berge, C. Principles of Combinatorics. Math, in Sc. and Eng. 72, Academic Press, New York, 1971. 
[3] Carlitz, L. Fibonacci Notes, 1. Zero-one sequences and Fibonacci numbers of higher order. Fibonacci 

Quart. 12 (1974), 1-10. 

[4] Fibonacci Notes, 5. Zero-one sequences again. Fibonacci Quart. (1977), 49-56. 

[5] Comptet, L. Analyse Combinatoire, vol. I. Presses Universitaires de France, Paris, 1970, pp. 32-34 and 

93-96. 

[6] Kaplansky, I. Solution of the "probleme des menages". Bull. Amer. Math. Soc. 49 (1953), 784-785. 
[7] Knuth, D. E. The Art of Computer Programming, vol. I: Fundamental Algorithms. Addison- Wesley, 
Reading, MA, 1968. 

[8] Montroll, E. W. Lattice Statistics, in Applied Combinatorial Mathematics. University of California 
Engineering and Physical Sciences Extension Series. E. F. Beckenbach editor. John Wiley, New York, 
1964. 

[9] Riordan, J. An Introduction to Combinatorial Analysis. John Wiley, New York, 1958. 
[10] Thompson, C. J. Mathematical Statistical Mechanics. Macmillan, New York, 1972. 

Department of Physics, University of Puerto Rico at Rio Piedras, PO BOX 23343, SAN 
JUAN PR 00931-3343 

E-mail address: tharrats@rrpac.upr.clu.edu 



